%% cal Q
%unit: m3/s
function Qa=Q_cal(Z,Qs,T,P)
P0=101.325*10^-3;
T0=293.15;
T=T+273.15;
P=P+P0;
Qa=Qs*Z*(P0/P)*(T/T0);

end